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ABSTRACT 


The  finite  element  code  Q3DFLO'81  was  evaluated  to  determine  its  suitability  for 
use  in  a  program  to  investigate  axial  compressor  tip  clearance  efiects.  The  code  was  first 
applied  to  Dring's  fully  comprehensive  benchmark  data  set  in  order  to  validate  the 
numerical  modelling  free  of  experimental  uncertainties.  It  was  then  applied  to  the  Naval 
Postgraduade  School  axial  research  compressor  with  which  tip  clearance  effects  were  to 
be  investigated  experimentally.  The  evaluation  identified  both  limitations  in  the  code 
and  limitations  in  the  data  which  could  be  obtained  in  the  experiment.  It  was 
reconuiiendcu  lliai  piovision  be  made  to  accommodate  peripherally  non-uniform  flow 
elTects  (blockage)  in  the  throughflow  code  calculation,  and  that  provision  be  made  to 
obtain  adequate  peripheral  fiow  surveys  in  the  experiment. 
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I.  INTRODUCTION 


The  efTect  of  finite  rotor  tip  clearance  on  the  operating  characteristics  of 
compressors  is  not  well  understood.  An  investigation  is  currently  underway  at  the 
Turbopropulsion  Laboratoiy  of  the  Naval  Postgraduate  School  (NPS)  to  explore  the 
interaction  of  the  tip  leakage  flow  with  the  rotor  passage  flow,  and  to  identify 
parameters  which  control  the  effect  of  tip  clearance  on  compressor  performance. 
Experiments  are  conducted  on  a  36  inch  diameter  low  speed  axial  flow  compressor  test 
rig.  The  compressor  has  three  repeating  stages  and  is  currently  configured  with  two 
stages  of  symmetrical  blading.  Measurements  of  the  compressor  performance  are 
obtained  using  fixed  probes  and  rakes,  a  flow  nozzle  and  a  torque  meter.  Plow  survey 
data  are  obtained  beween  blade  rows  using  pneumatic  probes.  Rotor  exit  conditions  are 
measured  using  hot-wire  probes.  Also,  using  high-response  semi-conductor  pressure 
transducers  and  conditional  sampling,  the  pressure  field  on  the  axisynunetric  surface 
across  the  rotor  tip  is  mapped.  .A  complete  description  of  the  three  dimensional  llowficld 
within  the  compressor  is  desired.  However,  for  each  type  of  flow  survey,  the  coverage 
of  the  compressor  annulus  is  limited  by  the  access  holes  provided  in  the  heasw 
compressor  case-wall.  Since  the  experimental  program  would  be  expected  to  develop  a 
modification  to  current  models  for  tip-clearance  effects  on  throughflow  development,  a 
computational  code  incorporating  a  case-wall  boundary  layer  calculation  was  needed 
with  which  to  obtain  predictions  of  the  experimental  conditions.  The  code  would  then 
provide  a  vehicle  in  which  to  incorporate  a  change  in  the  tip-flow  model. 

The  finite  element  code  Q3DFLO'81  was  examined  with  this  purpose  in  nund.  The 
code  was  installed  on  the  NPS  mainframe  IBM  370-3033  computer  in  19S3  by  its 
originator.  Professor  Charles  Hirsch,  and  was  used  sucessfully  at  that  time. 

The  purpose  of  the  present  evaluation  was  specifically  to  determine  the  suitability 
of  the  code  for  the  tip-clearance  application,  for  which  the  case-wall  boundary  layer 
modelling  would  be  a  significant  consideration.  The  evaluation  was  performed  in  two 
steps.  First,  the  code  was  applied  to  a  fully  documented,  'benchmark',  test  case  wherein 
completely  detailed  flow  field  information  was  available  in  addition  to  the  geometry  and 
the  controlled  boundary  conditions.  The  findings  from  this  effort  are  described  in 
Chapter  III.  Then,  with  a  developed  appreciation  of  the  degree  to  which  the  code 
described  actual  compressor  flow  conditions,  the  code  was  applied  to  the  NPS 
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compressor  geometry.  A  comparison  was  made  of  code  predictions  with  the 
pressure-rise  performance  measured  experimentally.  The  results  are  described  in 
Chapter  IV.  In  order  to  provide  a  comparison  with  measured  rotor- tip  pressure 
contours,  blade-to-blade  calculations  were  also  made  at  two  compressor  throughflow 
rates.  An  example  is  described  in  Chapter  V.  Conclusions  and  recommendations 
following  both  phases  of  the  evaluation  are  given  in  Chapter  VI.  In  general,  it  was 
found  that  the  most  significant  limitation  of  the  code  for  the  desired  application  lies  in 
the  absence  of  a  mechanism  to  input  peripherally  non-uniform  flow  effects,  which  are 
certainly  always  present  in  the  machine.  The  need  to  obtain  complete  peripheral  flow 
surveys  in  the  experimental  program  was  also  identified. 

Before  detailing  the  results  of  the  evaluation,  a  review  of  the  theoretical  background 
and  modeling  on  which  the  code  is  based,  is  given  in  Chapter  11. 
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II.  PRINCIPLES  OF  TURBOMACHINERY  THROUGHFLOVV 
MODELLING  AND  DESCRIPTION  OF  THE  CODE 


A.  BACKGROUND 

Q3DFLO'81  is  a  Finite  Element  code  developed  by  Prof  Charles  Hirsch  and  his 
coworkers  at  Vrije  University  in  Brussels  for  computing  the  complex  flowfield  in  a 
turbomachine  and  its  associated  axisymmetric  ducting.  The  formulation  of  the 
computational  model  of  the  flowfield  to  its  flnal  form  for  numerical  coding  is  described 
in  references  1  and  2.  Reference  1  describes  all  available  approaches  to  the  treatment 
of  turbomachinery  flowflelds  and  derisative  computational  methods,  including  those 
used  in  the  development  of  Q3DFL0'81.  Reference  2  is  a  summary  of  the  key  equations 
programmed  in  the  Q3DFLO'Sl  computational  code,  which  was  the  focus  of  the  present 
study. 

Since  the  two  references  were  intended  for  specialists  in  the  field  of  turbomachinery 
flow  modelling,  no  detailed  algebraic  operations  were  provided  from  one  step  to  another. 
For  a  newcomer  like  the  author,  it  was  not  easy  to  follow  the  logic  and  the  mathematical 
language  used  in  the  presentation.  To  complicate  matters,  the  equation  of  motion  can 
be  written  in  several  different  forms,  although  all  are  equivalent.  Typically,  for  a  flow 
in  an  annulus  region,  the  equation  is  normally  expressed  with  respect  to  an  absolute 
frame  of  reference  and  within  the  blade  to  blade  region  with  respect  to  a  relative  frame 
of  reference.  However,  in  the  Q3DFLO'81  formulation,  the  equation  of  motion  contains 
variables  from  both  the  absolute  frame  and  the  relative  frame,  in  the  same  equation! 

The  objective  of  this  chapter  therefore  is  to  review  the  development  of  the  equations 
which  are  the  basis  of  flowfield  computations  in  turbomachinery.  The  relationship 
between  the  various  forms  of  the  equation  of  motion  is  shown  and  the  detailed  algebraic 
operations  that  were  omitted  in  summarizing  the  development  of  the  Q3DFL0'81  in 
reference  2  are  described  here.  Detailed  derivation  of  equations  which  are  available  in 
texts  and  references  will  not  be  repeated.  In  such  situations,  only  the  key  equations  and 
their  physical  interpretations  are  presented.  It  is  hoped  that  this  will  give  similar 
newcomers  a  basic  understanding  before  attempting  the  advanced  materials  in 
reference  1. 
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B.  NON-AXISYMMETRIC  NATURE  OF  THE  REAL  FLOW 

Figure  1  shows  a  blade  row  with  radial  leading  and  trailing  edges  and  4  families  of 
surfaces;  namely,  surfaces  of  revolution  S„  ,  blade  surfaces  and  stream  surfaces  S,  and 
Si .  S„  is  the  axisymmetric  surface  generated  by  a  generatrix  m,  as  it  rotates  360’  about 
the  axis.  S,  is  simply  the  blade  surface.  Suppose  a  row  of  dyed  particles  is  initially  held 
along  the  edge  .A-P-C  and  another  row  along  5i-P-S,  at  the  inlet.  When  these  particles 
are  released,  it  is  found  that: 

•  The  row  of  particles  initially  along  A-P-C  does  not  follow  the  blade  surface  S^. 
They  instead  follow  the  stream  surface  Sj  contrary'  to  one's  expectation. 

•  Similarly,  the  row  of  particles  initially  along  does  not  follow  the 

axisymmetric  surface  S„.  Instead  they  follow'  the  stream  surface  S,. 

These  observations  have  been  verified  by  experiments.  A  flow  through  staiionaiy  or 
rotating  cascades  cannot  be  axisymmetric  if  the  flow  is  to  exert  a  moment  on  the  blades. 
[Ref.  3:  pp.  216] 

In  the  annulus  region  between  the  blade  rows,  the  flow  can  be  considered  to  be 
axisymmetric  after  the  wake  has  mixed  out.  In  the  formulation  of  through  How  analysis, 
the  flow  within  the  blade  row'  is  also  very  often  assumed  to  be  axisymmetric  to  simplify 
the  mathematical  model  i.e.  surfaces  S„  will  coincide  with  S,  and  with  S..  This 
assumption  is  also  adopted  in  the  following  derivation.  However,  this  assumption 
causes  discrepancies  in  the  velocity  components  and  flow  angles  between  experimental 
measurement  and  throughflow'  modelling  results.  These  discrepancies  can  be  corrected 
by  incorporating  a  factor  known  as  the  tangential  blockage.  This  is  dealt  with  in  the 
ne.xt  chapter. 

C.  EQUATION  OF  MOTION  IN  THE  ABSOLUTE  FRAME 

The  fundamental  equation  of  fluid  motion  can  be  expressed  as 

V,VV  =  -^-S7(gZ)  +  Fj-  {1} 

for  viscous,  compressible  or  incompressible  fluids. 

In  turbomachineiy  applications,  the  equation  of  motion  is  very  often  expressed  in 
terms  of  enthalpy  instead  of  pressure.  This  is  done  in  2  steps.  First,  use  the  vector 
identity 

F.VF  =  V(-^)- Fx  V  X  F  (2} 
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Second,  introduce  the  relationship  between  two  thermodynamic  states  along  the  path 
of  a  fluid  particle 

V/2-7V.S  =  -^  {3} 

Using  Eq.  (2)  and  Eq.  (3),  Eq.  (1)  becomes 

+  V(  -  K  X  V  X  E  =  TVS  -  V;?  -  V{gZ)  +  Ff  {4} 

Re-arranging, 

f  E  E^  “  ^  - 

-f-  \{h  +  ~  +  gZ)==  TVS  +  E  X  V  X  E  +  T;-  {5} 

The  total  enthalpy  is  defined  as 

//  =  /»  +  ~  +  gZ  (6] 

Using  Eq.  (6),  Eq.  (5)  becomes 


+ V//=  TVS-f  Ex  Vx  E-l-f.  {7} 

The  detailed  derivation  of  Eqn  (2)  and  (3)  is  given  in  reference  3.  For  a  steady  and 
isentropic  flow,  Eq.  (7)  reduces  to 


V//=ExVxE  {8} 

If  VT/  =  0,  the  fluid  is  either  non-rotational  or  the  vorticity  is  parallel  to  the  particle 
velocity. 

D.  EQUATION  OF  MOTION  WITHIN  A  BLADE  ROW 

Figure  2  shows  the  outlet  flow  of  a  constant  speed  rotating  rotor  with  a  relative 
velocity  profile  represented  by  the  vectors  IE,  to  A  stationary  observer  attached  to 
the  stator  blade  row  looking  at  the  rotor  outlet  will  see  a  non-steady  velocity  profile 
represented  by  the  vectors  f',  to  E^  which  vary^  in  both  magnitude  and  direction.  It  is 
clear  that  Eq.  (7)  is  not  useful  for  describing  the  flow  within  the  rotor  region,  as  the  term 
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cannot  be  neglected  nor  determined.  Hence,  it  is  more  useful  to  adopt  a  rotating 
frame  of  reference  that  rotates  at  the  same  speed  as  the  rotor  for  describing  the  flow 
within  the  rotor  blade  row.  [Ref  3:  pp.  110,111] 

E.  EQUATION  OF  MOTION  IN  THE  RELATIVE  FRAME 

The  equation  of  motion  with  respect  to  the  rotating  frame  can  be  obtained  directly 
from  Eq.  (1)  by  replacing  the  absolute  velocity  by  the  relative  velocity  using 

{9} 

and  taking  all  changes  with  respect  to  the  relative  system  (i.e.  Vj,  in  place  of  V  and  -:r“ 

^  0/ 

in  place  of-r-).  The  left  hand  side  (LHS)  of  Eqn  (1)  can  be  expressed  as 
c  t 

+  r .  V  r  =  -I-  ( (E  +  oj  X  R)  +  ( ir  +  w  X  R) .  v(  ir  +  w  X  R) 

Cl  Cl 

{10} 

=  -4— +  n  *v„fl-  +  2tox  ir  -  V;>( ) 

Cl  ^  "2 

and  Eq.  (1)  becomes 

(r.v,if+2coxir-v^(-^)  =  — ^ — v,(gZ)  +  F^  {ii} 

The  detailed  derivation  of  Eqn  (10)  is  given  in  reference  3  . 

Again  using  the  vector  identity  of  Eqn  (2)  and  the  thermodynamic  relationship  of 
Eqn  (3)  in  the  relative  frame.  Eqn  (11)  becomes 


CniV  -  (o^R^ 

-^  +  Vj,(-~)-  IVxV^x  iy+2coxJV-Vj,(^^) 

C I  ^  ^ 


=  T^rS  —  VrIi  —  Vjij(^Z)  +  /y 


Re-arranging, 


BrIV  rp2  (o^R^  - 

-|p  +  V;e(/i  +  -^4-gZ--^^^)=7V^5+  IFx(V;fX  W  +  2io)  +  Fj-  {13} 


With  the  rothalpy  defined  as 


{14} 


r  j/2  ^  n  ^ 

r  /  .  W  .  ^  CO  R 

l  =  h  +  —p  +  gZ-  — 2 — 


the  equation  of  motion  in  the  relative  frame  is  then 


c^JV 


Ip  +  V«/  =  TVftS  +  IV  X  (V«  X  W  +  2(o)  +  Ff 


{15} 


For  a  steady  and  isentropic  flow  within  the  rotor  and  with  the  rothalpy  constant,  Eq. 
(15)  reduces  to 


irx(V^x  lV+2^)  =  0  {16} 

Physically,  this  means  that  the  flow  is  rotational  with  respect  to  the  rotating  frame  and 
the  relative  vorticity  is  opposite  to  and  equal  to  twice  the  rotor  rotation. 

Flowever,  it  should  he  pointed  out  that  the  assumption  of  steady  relative  flow, 
though  frequently  made  in  turbomachinery  calculations,  is  true  only  if  the  inlet  flow  to 
the  rotor  is  uniform  in  the  tangential  direction;  i.e.  the  wakes  of  any  upstream  blade  rows 
must  have  mixed  out.  Hirsch  [Ref  1:  pp.  53]  states  that  in  a  real  flow,  "the  viscous 
interactions  and  unsteadiness,  which  are  always  present  to  some  degree  in 
turbomachiner.'  flow,  introduce  variations  in  rothalpy  along  the  flow  path." 

F.  EQUATION  OF  MOTION  USED  IN  Q3DFLO'81 

.An  alternative  form  of  the  equation  of  motion  for  steady,  relative,  three-dimensional 
flow  which  is  used  in  Q3DFL0'81  involves  both  relative  and  absolute  velocity 
components;  namely, 


V,i/ =  TV/fS  +  IV  X  (V  X  V)  +  Fj-  {17} 

Eq.  (17)  is  the  same  as  Eq.  (15)  for  a  steady  flow,  except  that  the  term  (V;;  x  IF  4-  2(1)) 
is  e.xpressed  in  terms  of  (V  x  F). 

G.  SIMPLIFIED  THROUGHFLOW  MODEL  OF  Q3DFLO'8I 

Figure  3  shows  the  characteristics  of  an  actual  flow'  within  the  blade  passage  of  a 
rotor.  The  flow  is  highly  three-dimensional  with  strong  comer  vortices  generated  at  the 
corners  between  the  blade  and  hub  and  betw'een  the  blade  and  shroud.  There  is  also  a 
secondary  flow'  surrounding  the  primary  core  flow'.  However,  the  Q3DFLO  S1 
throughflow  uses  a  simplified  axisymmetric  flow'  model,  and  the  three-dimensional 


7 


details  are  removed  (accounted  for)  by  taking  the  density  average  of  the  flow  over  the 
blade  spacing  at  each  spanwise  position. 

In  assuming  the  flow  to  be  axisymmetric  through  blade  row,  the  effects  of  the  blades 
themselves  must  be  introduced  into  the  equation  of  motion  as  a  body  force.  The 
equation  of  motion  (Eq.  17)  becomes 

-  JTx(Vx  K)- 7V/jS-V^/  +  fi,  +  ^  {18} 


H.  BASIC  FORMULATION  OF  Q3DFLO'81 

The  basic  steps  of  the  formulation  of  Q3DFLO'81  throughflow  model  are  described 
in  reference  2  .  The  throughflow  equation  known  as  the  Radial  Equilibrium  Equation 
that  is  solved  in  the  Q3DFLO'81  code  for  the  flow  within  the  rotor  is 


,  1  f"'/'  .  .  r  ,  1 

; — :: —  )  t-  (  —  ~ 

cr  pro  cr  oz  pro  cz 


)  = 


( ]V,  +  ir,  tan  A)(  4^  -  7-4f )  +  ( tan  r,  -  IV^)(  ^  -  T-|f ) 

C  f  C  r  C 


(19) 


If 


-2 


—  ( tan  — — - h  tan  j;  — — — ) 


cr 


cz 


Eq.  (19)  is  derived  from  Eq.  (18).  Details  of  the  derivation  are  given  in  Appendix  (.A). 

For  the  flow  within  the  stator,  Eq.  (19)  is  still  applicable.  As  the  stator  is  stationary, 
the  relative  velocity  components  are  replaced  by  the  absolute  velocity  components  and 
the  rothalpy  is  replaced  by  the  enthalpy. 

However,  for  the  flow  within  the  annular  region  without  blades  such  as  the  annular 
space  between  a  rotor  and  stator  blade  rows,  Eq.  (19)  is  not  applicable.  The  inviscid 
equation  of  motion  is  used  in  such  regions.  The  formulation  of  the  form  of  equation 
of  motion  that  is  solved  in  Q3DFL0'81  for  the  flow'  in  blade-free  annular  regions  is 
given  in  reference  2. 

Figure  4  shows  the  boundary  conditions  applied  to  the  throughflow  calculation. 
The  last  station  must  be  be  normal  to  the  meriodional  velocity. 

The  radial  equilibrium  equations  for  the  flow  within  the  blade  row  and  the  annular 
region  are  non-linear.  They  are  solved  by  the  finite  element  method  using  a  iterative 
procedure.  The  details  of  the  numerical  method  are  also  described  in  reference  2. 
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III.  ASSESSMENT  OF  THE  CODE  USING  A  BENCHMARK  DATA  SET 


A.  BACKGROUND 

The  set  of  radial  equilibrium  equations,  which  is  solved  in  Q3DFLO'81  is  only  one 
of  three  major  ingredients  in  the  code.  Q3DFLO'81  incorporates  an  endwall  boundary 
layer  calculation  which  (optionally)  can  be  coupled  into  the  radial  equilibrium  solution 
procedure.  The  flow  is  then  solved  iteratively  to  account  for  the  effect  of  boundary  layer 
and  tip  clearance  losses.  Like  all  other  throughflow  codes,  Q3DFLO'81  also  makes  use 
of  correlations  of  two  dimensional  cascade  data  for  the  loss  coefficient  and  flow  turning 
angle.  Consequently,  the  ability  to  predict  the  performance  of  an  axial  compressor  using 
Q3DFLO'Sl  depends  not  only  on  the  simplified  modelling  of  the  physical  flow  and  the 
accuracy  of  the  numerical  method  but  also  on  the  applicability  and  accuracy  of  the  two 
dimensional  cascade  data. 

Comparisons  of  predictions  using  Q3DFLO'81  with  experimental  data  on  axial  flow 
compressors  were  carried  out  at  the  Naval  Postgraduate  School  during  1983  and  1984. 
It  was  concluded  at  that  time  that  the  code  was  able  to  predict  the  general  behaviour 
of  the  flowfield  with  certain  discrepancies  [Ref  4].  However,  the  experimental  data  in 
each  application  were  not  complete  enough  to  enable  the  evalution  of  each  aspect  of  the 
code  independently,  such  as  the  axisymmetry  assumption  in  modelling  the  physical  flow 
and  the  use  of  two  dimensional  cascade  data  for  correlating  the  losses  across  the  blade 
row.  Furthermore,  the  spanwise  static  pressure  distribution,  wh’ch  is  the  most  sensitive 
indicator  of  accuracy  in  applying  throughflow  codes,  was  not  available  from  the 
experiments. 

In  1985,  under  the  sponsorship  of  the  Naval  Air  System  Command,  completely 
detailed  flowfield  surveys  obtained  in  a  five  feet  diameter  two  stage  axial  flow 
compressor  at  United  Technologies  Research  Center  (UTRC)  were  used  to  assess  the 
accuracy  and  utility  of  throughflow  codes  for  predicting  compressor  flow  fields  and 
compressor  performance.  The  relevant  parameters  of  the  flowfield  required  for  a 
comparison  with  code  predictions,  such  as  the  spanwise  variation  of  the  total  pressure 
loss  coefficient,  flow  turning  angle  and  tangential  aerodynamic  blockage,  were  obtained 
from  the  surveys.  Dring  and  Joslyn  [Ref  5:  pp.  56-64]  used  the  reduced  experimental 
data  in  a  throughflow  finite  element  code  developed  by  Habashi  and  Youngson.  They 
concluded  that  an  accurate  description  of  the  tangential  aerodynamic  blockage,  resulting 
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from  the  non-axisymmetr>'  of  the  physical  flow  and  the  presence  of  wakes,  was  essential 
if  the  flow  field  was  to  be  predicted  accurately.  The  detailed  set  of  measurements  was 
subsequently  made  available  as  a  benchmark  data  set  for  the  assessment  of  other 
throughflow  codes  in  reference  5. 

With  the  complete  survey  data  available  as  inputs,  the  Q3DFL0'81  code  can  be  run 
on  the  benchmark  compressor  without  accessing  the  two  dimensional  cascade  loss  and 
deviation  correlations  and  the  endwall  boundary  layer  subroutines.  This  enables  the 
effects  of  blade  wakes  and  accuracy  of  the  axisymmetric  flow  assumption  to  be  assessed 
by  running  the  code  with  and  without  accounting  for  the  tangential  aerodynamic 
blockage.  In  the  process,  the  accuracy  of  the  code's  numerical  method  can  be  exantined. 

The  benchmark  axial  flow  compressor  had  extensive  hub  corner  stall  in  the  second 
stage  stator.  Although  cascade  data  are  not  appropriate  for  correlating  the  losses  and 
exit  deviation  angle  across  the  blade  row  of  compressor  with  extensive  region  of  stall,  it 
is  interesting  to  see  how  diflerent,  quantitively,  the  predictions  will  be  using  the 
correlation  subroutine  in  Q3DFLO'8l.  which  is  based  on  two  dimensional  cascade  data, 
compared  to  the  use  of  actual  survey  data.  This  interest  arises  because  Dring  has 
observed  the  occurrence  of  hub  corner  stall  in  all  the  configurations  tested  at  UTRC. 

B.  BLOCKAGE  FACTORS  IN  THE  Q3DFLO'81  FORMULATION 

There  arc  2  types  of  blockage  involved  in  throughflow  modelling;  namely,  endwall 
blockage  and  tangential  (two  dimensional)  blockage  [Refs.  2,6].  Endwall  blockage 
accounts  for  the  reduction  in  throughflow  area  due  to  the  displacement  thickness  of  the 
boundary  layer  along  the  shroud  and  hub  surfaces.  Q3DFL0'81  accounts  for  the  effect 
of  endwall  blockage  by  applying  the  endwall  boundary  conditions 

4>  =  - 
and 

at  the  hub  and  case  walls  respectively.  Here,  m  is  the  massflow  rate  and  5'^  and  (52  are 
the  boundary  layer  displacement  thicknesses  at  the  hub  and  shroud  respectively. 

The  tangential  blockage,  (1  -  b),  accounts  for  the  reduction  in  the  effective  annulus 
area  in  the  circumferential  direction.  The  tangential  block •’ge  varies  as  a  function  of 
radius.  Dring  [Ref.  5:  pp.  1 1]  states  that  "the  blockage  (1  -  b)  reflects  the  fraction  of  the 
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circumference  not  available  to  the  inviscid  throughflovi'."  Tangential  blockage  can  be 
mechanical  or  aerodynamic.  Mechanical  blockage  occurs  within  the  blade  row  due  to 
the  airfoil  tangential  thickness.  Aerodynamic  blockage  occurs  both  within  the  blade  row 
as  well  as  in  the  annulus  region  between  blade  rows.  Within  the  blade  row.  the 
aerodynamic  blockage  results  from  the  displacement  thicknesses  of  the  boundary  layers 
along  the  blade  surfaces  and  the  departure  from  axisymmetry.  In  the  annulus  region, 
the  aerodynamic  blockage  is  due  to  blade  wakes  and  non-axisymmetries  such  as  corner 
stall. 

Q3DFLO'81  accounts  for  the  tangential  blockage  by  defining  the  a.xial  (IT,)  and 
radial  ( H ',)  velocity  components  in  its  governing  equation  as 


and 


I  c'A 
prb  cz 


However.  Q3DFL0'81  only  computes  the  tangential  blockage  (1  -  b)  due  to  the 
mechanical  blockage  within  the  blade  rows,  and  does  not  include  aerodynamic  blockage. 
There  is  no  provision  for  introducing  the  aerodynamic  blockage  as  input  data. 

C.  ASSESSMENT  METHODOLOGY 

.Measured  values  of  deviation  angle  and  loss  coefficients  can  be  used  as  inputs  to  the 
code.  Consequently,  the  applicability  of  the  correlation  expressions  for  the  losses  and 
deviation  angle  that  are  contained  in  Q3DFLO'Sl.  to  the  prediction  of  the  performance 
of  a  compressor  with  a  large  stall  region,  can  be  assessed.  First,  the  code  can  be  run 
inputting  measured  quantities  from  the  benchmark  data  set.  and  then  run  again  using 
the  correlation  subroutines.  Both  runs  will  use  the  same  geometrical  data  and  will  have 
the  same  inlet  conditions. 

An  assessment  of  the  eflcct  of  including  tangential  aerodynamic  blockage  cannot 
be  done  in  a  straight  forward  manner  since  Q3DFLO'81  does  not  have  provisions  for 
introducing  tangential  aerodynamic  blockage  as  input  data.  What  can  be  done  is  to  use 
an  average  value  for  the  tangential  aerodynamic  blockage  (given  by  reference  5  on  page 
21  as  a  mass-average  value)  and  equate  it  to  the  equivalent  endwall  blockage,  which  the 
code  can  accept  as  an  input. 
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This  approach,  although  not  physically  correct,  does  give  an  indication  of  the  effect 
of  the  tangential  aerodynamic  blockage  on  the  throughllow  prediction.  In  a  broad 
sense,  both  types  of  blockages  are  similar  in  their  effect.  They  both  increase  the 
meriodional  velocity  component  in  the  ihroughflow.  From  the  results  presented  by 
Dring  and  Joslyn  [Ref  5;  pp.  49,99,102],  one  can  conclude  that,  in  the  absence  of  a 
detailed  description  of  the  radial  variation  in  the  tangential  aerodynamic  blockage,  a 
uniform  tangential  aerodynamic  blockage  (an  average  value)  produces  much  better 
results  than  by  not  including  the  tangential  aerodynamic  blockage  at  all. 

Three  cases  were  run  on  the  benchmark  compressor  using  Q3DFLO’81.  Cases  1  and 
2  were  run  using  the  measured  values  of  loss  and  deviation  angle.  However,  Case  1  was 
run  with  no  tangential  aerodynamic  blockage  and  Case  2  with  an  equivalent  tangential 
aerodynamic  blockage  added  to  the  case-wall.  The  results  of  Cases  1  and  2  are 
compared  with  the  benchmark  compressor  experimental  results  to  assess  the  effect  of 
including  tangential  aerodynamic  blockage  on  the  prediction.  Case  3  was  run  using  the 
Q3DFL0'81  correlation  subroutines  and  with  the  same  equivalent  tangential 
aerodynamic  blockage  added  to  the  case-wall  as  for  Case  2.  The  Case  3  prediction  was 
compared  to  those  of  Cases  1  and  2.  In  all  the  three  cases,  the  endwall  boundaiy  layer 
computation  in  the  code  was  not  used. 

D,  INPUT  DATA 

The  annulus  of  the  benchmark  compressor  was  cylindrical.  The  radius  of  the  hub 
and  tip  were  24  ins  and  30  ins  respectively.  The  airfoils  of  the  rotor  and  stator  were 
NACA  65  series  with  circular  arc  camber  lines.  Figure  5  and  Figure  6  show  a  schematic 
and  the  computational  mesh  for  the  benchmark  compressor  respectively.  Onl}  the 
second  stage  rotor  and  stator  were  modelled.  This  was  because  the  detailed 
measurements  were  provided  for  the  second  stage  and  furthermore,  Q3DFLO’81  is 
limited  by  the  griJ  capacity  and  could  not  model  the  whole  compressor  and  maintain 
adequate  resolution  near  the  walls.  A  total  of  15  stations  was  used.  Station  1.  the  inlet 
station  to  the  second  stage,  is  located  mid-way  between  the  first  stage  stator  trailing  edge 
and  the  second  stage  rotor  leading  edge.  Stations  2  and  6  correspond  to  the  second  stage 
rotor  leading  and  trailing  edges  respectively,  and  Stations  8  and  12  to  the  second  stage 
stator  leading  and  trailing  edges  respectively.  T6e  code  predictions  and  benchmark 
measurements  were  compared  at  Stations  7  and  13,  which  correspond  to  Stations  4  and 
5  in  reference  5. 
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A  total  of  1 1  axial  grid  lines  (which  is  the  maximum  permitted  by  Q3DFL0'81)  was 
used.  Ihe  streamlines  were  also  clustered  towards  the  shroud  by  using  a  repartition 
factor  of  0.1,  This  emphasis  on  the  tip  region  was  necessitated  by  the  abrupt  variation 
in  the  benchmark  loss  coefficient  near  the  casewall.  Figure  7  shows  the  distributions 
of  the  velocities  and  total  pressure  as  a  function  of  radius  at  Station  1,  which  were  input 
for  all  three  cases.  Figure  8  and  Figure  9  show'  the  distributions  of  the  benchmark  loss 
coefficient  and  exit  flow  angle  as  a  function  of  radius  for  the  second  stage  rotor  and 
stator,  which  were  input  for  Cases  1  and  2. 

Several  points  are  noted.  The  benchmark  data  are  given  in  terms  of 
non-dimensional  quantities.  For  use  in  Q3DFLO'81,  the  data  were  converted  to  the 
dimensional  form  in  metric  units.  A  Standard  Day  of  15°  C  and  1.01325  bar  was  chosen 
as  the  inlet  condition  of  tlie  compressor  to  calculate  the  mass  flow  rate  and  the 
velocities.  The  loss  coefficients  given  in  the  benchmark  data  are  specified  with  respect 
to  the  mid-span  wheel  speed.  For  Q3DFLO'81,  the  loss  coeflicients  were  re-calculated 
with  respect  to  the  benchmark  rotor  inlet  relative  velocity  and  the  stator  inlet  absolute 
velocity  at  their  specified  radial  positions,  for  the  rotor  and  stator  respectively.  In 
Q3DFLO'Sl,  the  aerodynamic  input  data  are  limited  to  10  radial  points  for  each 
variable.  Since  IS  to  20  points  were  specified  in  the  benchmark  data  set  for  each 
variable,  ten  points  were  carefully  selected  such  that  they  represented  the  macro 
variation  as  closely  as  possible. 

E.  EFFECT  OF  TANGENTIAL  AERODYNAMIC  BLOCKAGE  ON  THE 
THROLGHFLOW  PREDICTION 

Figure  10  through  Figure  17  compare  the  Q3DFL0'81  code  predictions  for  Cases 
1  and  2  with  the  experimental  results.  The  comparisons  are  made  at  Stations  7  and  13 
for  the  following  parameters: 

•  Total  pressure  coefficient 

•  Relative  total  pressure  coefficient 

•  Static  pressure  coefficient 

•  Axial  velocity  component 

•  Absolute  flow  velocity 

•  Relative  flow  velocity 

•  Relative  flow  yaw  angle 

•  .Absolute  flow  yaw  angle 
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Q3DFLO'81  outputs  the  pressure  and  flow  velocity  in  metric  units;  namely,  bar,  and 
meter  per  second.  These  parameters,  shown  plotted  in  Figure  10  to  Figure  17,  were 
non-dimensionalised  using  the  mid-span  wheel  speed. 

At  Station  7,  Figure  10  through  Figure  17  show  that  the  predictions  of  Cases  '  and 
2  are  very  similar.  It  can  be  seen  that  Case  2  agrees  slightly  better  than  Case  1  with  the 
benchmark  data  in  both  magnitudes  and  overall  trends. 

At  Station  13,  Figure  12,  Figure  13  and  Figure  16,  however  show  that  there  are 
significant  differences  in  the  predictions  of  static  pressure,  axial  velocity  and  relative  flow 
yaw  angle  between  Case  1  and  Case  2.  The  prediction  of  Case  2  agrees  very  well  with 
the  benchmark  data,  but  not  Case  1.  Case  1  shows  appreciably  higher  static  pressure 
and  relative  flow  yaw  angle  and  lower  axial  velocity. 

The  agreement  at  Station  7  and  disagreement  at  Station  13  between  Case  1  and  Case 
2  are,  almost  certainly,  the  result  of  the  tangential  aerodynamic  blockage  introduced  at 
these  stations  in  Case  2.  The  tangential  aerodynamic  blockage  is  1.4  percent  at  Station 
7  compared  to  4.7  percent  at  Station  13.  For  a  constant  massflow  rate,  the  absence  of 
1,4  percent  area  blockage  at  station  7  in  Case  1  causes  a  slight  reduction  in  the  axial 
velocity  which  results  in  only  a  slight  increase  in  static  pressure  and  absolute  flow  yaw 
angle  as  can  be  seen  in  Figure  12  and  Figure  17.  However,  the  effect  of  ignoring  the 
4,7  percent  blockage  at  Station  13  in  Case  1  is  clearly  significant.  It  causes  a  marked 
reduction  in  the  predicted  axial  velocity,  which  results  in  a  significant  increase  in  the 
static  pressure  and  relative  flow  yaw  angle. 

F.  CODE  PREDICTION  USING  TVVO  DIMENSIONAL  CASCADE  DATA 

Figure  10  through  Figure  17  also  show  the  predictions  for  Case  3.  which  uses  the 
code's  correlations  of  cascade  data,  rather  than  measured  data,  for  computing  loss  and 
deviation  angle.  The  Q3DFL0'81  endwall  boundary  layer  part  of  the  code  was  not  used 
because  of  the  limited  resolution  in  the  benchmark  data  and  difficulties  in  starting  the 
endwall  boundary  layer  calculation  in  the  middle  of  the  compressor.  Clearly,  the  most 
meaningful  evaluation  of  the  accuracy  of  the  cascade  correlation  will  come  from  a 
comparison  of  measurements  and  predictions  in  the  core  flow,  say  from  25  percent  to 
75  percent  span  from  the  hub. 

Within  the  core  flow,  the  Case  3  predictions  of  the  velocities  and  flow  angles  at  both 
Stations  7  and  13  generally  follow'  the  trends  in  the  benchmark  experimental  results. 
However,  they  are  not  as  good  as  the  predictions  obtained  in  Cases  1  and  2.  The 
maximum  disagreement  in  the  flow  angle  is  seen  to  be  less  than  three  degrees  and  in  the 
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flow  velocity,  less  than  ten  percent.  The  code  also  computed  much  higher  total  and 
static  pressures,  especially  at  Station  13.  This  would  be  expected  perhaps,  because  of  the 
higher  losses  resulting  from  the  presence  of  hub  corner  stall. 

G.  OBSERVATIONS 

Comparison  of  the  Q3DFLO'81  predictions  for  Cases  1  and  2  clearly  demonstrate 
the  need  to  include  the  tangential  aerodynamic  blockage  resulting  from  the  wakes,  or 
an  equivalent  measure  of  the  departure  from  axisymmetry,  in  computing  the  compressor 
throughflow.  The  good  agreement  of  the  predictions  of  Case  2  with  the  benchmark 
experimental  results,  is  considered  to  validate  the  code's  numerical  method. 

Results  from  Case  3  also  show  that  it  is  necessary  to  know  in  advance  whether  the 
compressor  being  modelled  has  an  extensive  stalled  region  of  flow.  If  the  compressor 
has  an  extensive  region  which  is  stalled,  the  subroutines  contained  in  Q3DFLO'81  for 
calculating  the  total  pressure  losses  and  deviation  angles  will  not  give  accurate  results. 

However,  in  general,  given  the  correct  inputs,  the  simplified  throughflow  model  is 
able  to  predict  correctly,  both  the  levels  and  the  spanwise  distributions  of  the  flow  at 
each  station  in  an  axial  flow  compressor. 
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IV.  APPLICATION  TO  A  MULTI-STAGE  COMPRESSOR 

A.  BACKGROUND 

The  three-stage,  axial  flow  compressor  test  rig  at  the  Turbopropulsion  Laboratory 
of  the  Naval  Postgraduate  School  (NPS)  is  presently  configured  for  an  experimental 
program  to  determine  the  effect  of  rotor  tip  clearance  on  the  efficiency  and  flow  behavior 
in  an  axial  flow  machine.  Figure  18  shows  a  radial  section  of  the  compressor.  It  has 
three  repeating  stages  of  symmetrical  blading.  It  also  has  a  row  of  inlet  guide  vanes  to 
provide  the  required  pre-rotation  to  the  first  stage  rotor.  The  stator  provides  the  same 
pre-rotation  to  each  following  stage.  The  design,  geometry  and  construction  of  the 
compressor  are  described  in  references  7  and  8. 

Q3DFL0'81  was  applied  to  analyse  the  machine  performance  and  flow  through  the 
test  rig.  The  purpose  was  to  see  whether  Q3DFLO'8I  could  be  used  to  support  the 
experimental  research  activities  on  the  effect  of  rotor  tip  clearance.  Furthermore,  if 
Q3DFL0'81  was  shown  to  predict  the  performance  of  the  test  rig  correctly,  then  the 
predicted  throughflow  solution  could  be  used  as  input  for  the  blade-to-blade  analysis, 
which  is  also  available  in  Q3DFLO'81. 

The  present  chapter  deals  with  the  application  of  the  throughflow  code  contained 
in  Q3DFL0’81  to  the  NPS  compressor  test  rig.  In  the  process,  some  practical  aspects 
of  throughflow  modelling  using  Q3DFLO'81  were  uncovered.  The  configuration  of  the 
test  rig  modelled  here  consisted  of  the  IGV's  and  only  the  second  and  third  stages,  since 
this  was  the  configuration  being  used  in  the  experimental  investigation  of  the  effects  of 
tip  clearance  change. 

B,  EXPERIMENTAL  SET-UP 

Figure  18  also  shows  the  experimental  arrangement  that  was  used  to  measure  the 
overall  stage  pressure  rise  of  the  compressor.  The  overall  static  pressure  rise  was 
measured  using  two  static  pressure  taps  located  on  the  shroud  1.6  inches  ahead  of  the 
second  stage  rotor  and  two  taps  at  the  exit  of  the  third  stage  stator.  Each  pair  of  taps 
were  spaced  about  90°  apart. 

The  mass  flow  rate  was  calculated  from  an  inlet  bellmouth  measurement,  involving 
four  impact  pressure  probes  and  four  static  pressure  taps  spaced  equally  apart,  and  one 
thermocouple  sensor. 
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For  the  purpose  of  prescribing  the  inlet  conditions  for  the  code,  the  flow  at  the  inlet 
to  the  IGV  was  measured.  The  total  temperature  was  assumed  to  be  uniform  and  the 
same  as  for  the  bellmouth  measurement.  The  spanvvise  distribution  of  the  total  pressure 
was  measured  using  a  nine  probe  total  pressure  rake  located  0.87  inches  downstream  of 
the  exit  plane  of  the  strut  which  was  1.75  inches  upstream  of  the  IGV.  The  nine  probes 
were  aligned  along  a  radial  line  positioned  almost  mid-way  between  two  adjacent  struts. 
The  two  probes  nearest  the  endwalls  were  0.2  inches  and  0.25  inches  away  from  the  hub 
and  shroud  respectively.  The  static  pressure  was  measured  at  the  shroud  and  the  hub 
using  two  static  taps,  and  at  mid-channel  using  a  pilot  static  pressure  probe.  The 
velocity  was  assumed  to  be  axial  and  its  magnitude  was  calculated  from  the  dynamic 
pressure  distribution. 

All  measurements  were  recorded  using  a  Hewlett  Packard  3053  Data  Acquisition 
System.  A  description  of  the  data  acquisition  system,  the  instrumentation  and  the 
measurement  uncertainties  are  given  in  reference  9.  The  uncertainty  of  pressure  rise 
coelTicient  and  flow  coefficient  was  established  from  repeated  measurements  to  be  0.5 
percent  and  0.6  percent  respectively. 

C.  EXPERIMENTAL  MEASUREMENTS 

Experimental  measurements  were  obtained  at  a  fixed  compressor  speed  for  four 
mass  flow  rates.  For  performance  mapping,  the  mass  flow  rate  is  expressed  in  terms  of 
the  dimensionless  flow  coefiicient  (0)  defined  as 


^  n-  A  r 

where  m  is  the  massflow  rate  determined  from  the  bellmouth  measurement,  and 
are  the  density  and  annulus  area  respectively  at  the  IGV  inlet,  and  U„  is  the  mid-span 
peripheral  speed  of  the  rotor. 

The  floAV  coefficients  so  obtained  were  0.61,  0.67,  0.70  and  0.75.  Table  1  tabulates 
the  pertinent  experimental  measurements  and  calculated  parameters.  The  static  pressure 
rise  (AP,)  across  the  rotors  and  stators  is  expressed  by  the  dimensionless  quantity,  FI,, 
which  is  defined  as 


n,= 


2  Pioy^ 
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The  variation  of  H,  with  0  is  the  compressor  'characteristic'.  In  the  following  section, 
the  measured  characteristic  is  compared  with  the  Q3DFL0'81  code  prediction. 


Table  I.  EXPERIMENTAL  FLOW  AND  PRESSURE  COEFFICIENTS 


4> 

0.612 

0.670 

0.701 

0.748 

n. 

1.050 

0.952 

0.887 

0.799 

T,.:rR) 

534.56 

533.55 

534.21 

534.88 

Pot  (inw) 

409.56 

409.51 

409.49 

409.46 

Pos  (inw) 

407.53 

407.06 

406.79 

406.37 

Pm,  (Ibm/F) 

0.074223 

0.074278 

0.074137 

0.073967 

P'„,,  (ft  sec) 

95.50 

104.88 

110.21 

118.03 

i!'-.,,;  (ins) 

36.0 

36.0 

36.0 

36.0 

Aum 

7.06S6 

7.0686 

7.0686 

7.0686 

m  (Ibm  sec) 

50.104 

55.066 

57.753 

61.712 

ZXR) 

533.56 

532.55 

533.21 

533.88 

P,  (  inw) 

392.91 

393.50 

393.83 

394.40 

(Ibm/)’) 

0.071748 

0.071992 

0.071963 

0.071978 

P„  (inw) 

392.58 

393.11 

393.35 

393.87 

P,.  (inw) 

407.00 

406.24 

405.59 

404.88 

AP,  (ijnv) 

14.42 

13.13 

12.24 

11.01 

X  (rpm) 

1609.6 

1609.6 

1610.6 

1610.8 

R„  (ft) 

1.20 

1.20 

1.20 

1.20 

(ft  sec) 

202.27 

202.27 

202.39 

202.42 

The  measured  spanwise  distributions  of  total  and  static  pressure  and  the  velocity 
calculated  from  tlieni,  are  shown  in  Figure  19.  The  static  pressure  is  observed  to  be 
skewed  linearly  from  hub  to  shroud  based  on  the  shroud,  mid-span  and  hub 
measurements.  The  (axial)  velocity  was  calculated  from  the  measurements  using 

_  V  igv  ~ 


where  the  density  (p,^,)  in  pounds  mass  per  cubic  foot  at  the  IGV  inlet  was  obtained 
using 


f'igy 


igv  518.69 
L.igy  ^  ^^08.2 


X  0.07647 
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D.  MODELLING  OF  THE  COMPRESSOR 

The  NPS  3-stage  axial  compressor  has  a  cylindrical  annulus.  The  radius  of  the  hub 
and  tip  are  10.8  inches  and  18  inches  respectively.  The  airfoils  of  the  rotor  and  stator 
are  both  of  the  C  series  type  with  a  modified  thickness  distribution.  The  compressor 
was  modelled  from  the  measurement  plane  upstream  of  the  IGV  to  the  station 
downstream  of  the  third  stage  stator.  Figure  20  shows  the  computational  mesh,  which 
had  a  total  of  25  stations  and  10  axial  grid  lines.  Grid  lines  were  clustered  towards  the 
shroud  by  specifying  a  repartition  factor  of  0.1.  As  the  total  number  of  nodes  then 
exceeded  the  maximum  of  600  nodes  allow'ed  in  the  Q3DFL0'81  code,  the  analysis  was 
done  in  two  parts.  The  first  part  analysed  the  flow  from  Station  1,  ahead  of  the  IGV, 
to  a  station  immediately  downstream  of  the  third  stage  rotor.  The  prediction  at  this 
station  was  then  used  as  input  to  a  second  calculation  from  the  inlet  of  the  third  stage 
stator  to  the  measurement  station  downstream  of  the  stator  exit. 

The  input  data  required  by  the  code  for  the  inlet  condition  at  Station  1  were  the 
mass  flow  rate,  the  spanwise  distributions  of  the  velocity  components,  the  total 
temperature  and  the  total  pressure. 

E.  EFFECT  OF  THE  STRUT  WAKES 

The  boundary  layer  displacement  thickness  at  the  shroud  was  measured  to  be  about 
0.13  inches  in  a  separate  study  by  Tarigan  [Ref.  10).  For  the  code  prediction,  a  nominal 
endwall  blockage  of  2.6  percent  was  specified  at  Station  1,  which  is  equivalent  (in  area) 
to  a  0.15  inches  displacement  thickness  at  the  shroud.  The  higher  blockage  factor  was 
to  account  for  the  thicker  boundary  layer  at  the  junction  of  the  strut  and  the  endwalls 
and  the  negligibly  thin  boundary  layer  at  the  hub. 

A  disagreement  with  the  mass  flow  measurement  obtained  from  the  inlet  bellmouth 
was  noted  when  the  spanwise  velocity  distribution  at  the  IGV  inlet  was  integrated  to 
determine  the  mass  flow'  rate.  The  data  in  Table  2  show  that  the  massflow  rates 
determined  from  the  measurements  at  Station  1  were  about  3.5  percent  higher  than  the 
bellmouth  measurement  for  all  four  flow  coefficients. 
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Table  2.  BELLMOUTH  VS  INLET  RAKE  MEASUREMENTS  OF  FLOW 
RATE 


Examination  revealed  that  the  measurement  plane  at  Station  1  was  about  21  percent 
of  chord  behind  the  strut  trailing  edge.  There  were  a  total  of  six  struts  equally  spaced. 
Figure  21  shows  the  profile  of  the  strut  which  had  a  chord  of  4.12  inches  and  maximum 
thickness  of  0.813  inches.  It  was  known  that  strut  wakes  were  present  at  the 
measurement  plane,  although  they  did  not  intercept  the  rake  sensors.  The  thickness  of 
the  wake  had  not  been  measured  since  the  mechanical  arrangement  of  the  casing  did  not 
permit  a  circumferential  survey  of  the  flow.  However,  two  static  pressure  taps  on  the 
shroud,  located  between  the  IGV  outlet  and  the  2nd  stage  rotor  inlet,  did  indicate  an 
increase  in  the  static  pressure,  as  the  flow  proceeded  downstream.  This  could  be 
explained  by  the  diffusion  of  the  strut  wakes  in  the  annulus.  The  presence  of  the  strut 
wakes  also  explains  the  apparent  increase  in  mass  flow  rate  indicated  by  the  rake  at 
Station  1.  The  strut  wakes,  were  effectively  a  two  dimensional  blockage,  causing  a 
reduction  in  the  circumferential  flow  area  and  an  increase  in  the  local  flow  velocity  from 
hub  to  shroud  at  the  IGV  measurement  plane. 

F.  ADAPTATION  OF  THE  INLET  FLOW  CONDITION 

The  displacement  thickness  of  the  strut  wakes  could  be  expected  to  be  fairly  uniform 
in  the  core  of  the  flow  but  not  at  the  corners  between  the  strut  and  the  endwalls.  For 
numerical  modelling,  it  was  reasonable  to  assume  that  it  was  uniform  from  hub  to 
shroud.  The  endwall  boundar>'  layer  displacement  thickness,  the  mass  flow  rate 
determined  from  the  bellmouth  measurement  and  the  spanwise  axial  velocity  measured 
at  Station  I  were  sufficient  to  calculate  a  unique  displacement  thickness  of  the  strut 
wake  which  gave  the  measured  mass  flow  rate. 

The  wake  displacement  thickness  can  be  viewed  as  tangential  blockage. 
Unfortunately,  the  Q3DFL0'81  code  was  unable  to  accept  a  tangential  blockage  factor 
as  input  data.  Consequently,  it  was  assumed  that  the  strut  wakes  were  completely 
diffused  at  constant  radii  before  entering  the  IGV,  and  this  diffusion  process  was 


calculated  externally  from  the  Q3DFLO'81  code.  Appendix  (B)  describes  the  diffusion 
calculation  for  the  strut  wakes.  The  spanwise  distributions  of  the  diffused  (or  'mixed 
out')  axial  velocity  and  total  pressure  were  then  used  as  boundary  conditions  for  the 
code,  at  Station  1. 

One  could  expect  that,  by  assuming  the  wakes  had  diffused  before  entering  the  IGV, 
instead  of  letting  them  be  swallowed  by  the  IGV's  and  then  diffused  gradually  to  have 
negligible  thickness  at  Station  8,  the  prediction  of  the  axial  velocity  component  would 
differ  from  the  actual  value  immediately  downstream  of  the  IGV.  But  at  Station  8,  and 
beyond  where  it  was  expected  that  the  strut  wakes  would  have  diffused  completely,  the 
code  prediction  would  be  valid.  A  comparison  of  the  code  predictions  with  the 
experimental  measurements  for  the  static  pressure  at  the  shroud  at  Station  8  would  give 
a  good  indication  of  the  validity  of  the  above  reasoning. 

G.  RESULTS  FOR  THE  ADAPTED  INLET  FLOW 

The  endwall  blockage  and  wake  thickness  of  each  strut  at  Station  1  required  for  the 
code  calculation  to  match  the  measured  static  pressure  at  Station  S,  are  shown  in 
Table  3. 


Table  3.  DISPLACEMENT  AND  WAKE  THICKNESSES  AT  STATION  1 


Flow  Coeficient 

0.612 

0.670 

0.701 

0.748 

Endwall  Blockage  ^  Stn  1  (''o) 

2.6 

2.6 

2.6 

2.6 

Wake  Thickness  (inch) 

IB3H 

0.528 

0,528 

0.413 

Computed  P,  a  Stn  8  (inw) 

392.66 

393.10 

393,4d 

393.80 

Measured  P.  a  Stn  S  (inw) 

392.58 

393.10 

393,35 

393.87 

It  was  found  that  a  wake  displacement  thickness  of  0.528  inches  and  an  endwall 
blockage  factor  of  0,970  produced  a  match  with  the  measured  static  pressure  at  Station 
8  fairly  well.  The  ratio  of  the  wake  displacement  thickness  to  the  strut  maximum 
thickness  was  0.59. 

Figure  22  shows  the  spanwise  distributions  of  total  and  static  pressures  and  axial 
velocity  after  the  diffusion.  It  can  be  seen  that  the  total  pressure  was  hardly  affected. 
As  expected,  there  was  an  appreciable  increase  in  the  static  pressure  of  about  0.5  inches 
of  w.-’ter  due  to  a  3.5  percent  reduction  in  velocity. 
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H.  COMPARISON  OF  THE  OVERALL  PERFORMANCE 

Q3DFLO'81  was  used  to  predict  the  overall  performance  at  each  of  the  four  flow 
coefficients  both  with  and  without  strut  wake  diffusion  at  Station  1.  The  computed 
variation  of  O,  vs  4>  is  shown  in  Figure  23  denoted  as  follows: 

•  Case  1:  Experimental  Measurement. 

•  Case  2:  Q3DFLO'8I  Prediction  with  Correction  for  Strut  Wake  Blockage. 

•  Case  3:  Q3DFL0'81  Prediction  without  Correction  for  Strut  Wake  Blockage. 

•  Case  4:  .Meanline  Analysis  (Ref  9) 

It  should  be  noted  that  the  flow  coefficient  of  Case  3  was  based  on  the  actual  mass  flow 
rate  rather  than  the  apparently  higher  mass  flow  rate  measured  at  the  IGV  inlet. 

Comparing  Cases  2  and  3,  Case  2  shows  significantly  higher  pressure  rise  coefficients 
than  does  Case  3.  This  is  probably  for  two  reasons.  First,  the  reduced  axial  velocity 
vvhich  results  from  wake  diffusion  causes  an  increase  in  the  blade  loading.  The  IGV^  exit 
flow  angle  calculated  by  the  code  is  independent  of  the  value  of  the  flow  coefficient.  The 
code  predicted  essentially  the  same  absolute  flow  angles  at  all  stations  from  the  exit  of 
the  IGV  at  Station  4  to  the  inlet  of  the  first  rotor  at  Station  9,  both  for  Cases  2  and  3. 
With  the  same  absolute  flow  angle,  and  lower  a.xial  velocity,  Case  2  involves  a 
correspondingly  lower  absolute  tangential  velocity  upstream  of  the  rotor.  This  would 
cause  an  increase  in  the  relative  tangential  velocity  component,  which  results  in  an 
increase  in  the  blade  loading  of  the  first  rotor.  Second,  the  total  pressure  loss  coefficient 
calculated  by  the  Q3DFLO'Sl  subroutine  for  Case  2  was  significantly  lower  than  for 
Case  3  for  all  blade  rows.  This  suggested  that  Case  2  involved  more  favorable  incidence 
angles  than  Case  3. 

Comparing  Case  2  with  Cases  1  and  4,  Case  2  showed  better  agreement  with  Case 
1  (measurements')  at  the  flow  coefficient  of  0.612  but  departed  by  an  increasing  margin 
for  flow  coefficient  above  0.67.  Since  the  design  flow  coefficient  was  about  0.64,  the 
prediction  was  better  near  the  design  condition  than  at  off-design  conditions,  which 
might  be  expected. 

I.  EFFECT  OF  THE  IGV  EXIT  FLOW  UNDERTURNING 

It  was  suspected  that  the  difference  in  the  pressure  rise  coefficient  between  the 
experimental  measurements  and  the  predictions  at  all  conditions  might  be  the  result  of 
underturning  of  the  IGV  exit  flow,  which  had  been  observed  in  the  compressor. 
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Waddell  (Ref.  11]  had  measured  the  IGV  exit  flow  angle  and  had  established  that 
it  was  independent  of  the  flow  coeflicient.  He  showed  that  the  turning  angle  was  smaller 
than  the  design  intent  over  a  large  part  of  the  span.  The  Q3DFLO'81  prediction  of  the 
IGV  exit  flow  angle  for  Case  2  is  compared  with  Waddell's  measurement  in  Figure  24. 
It  can  be  seen  that  the  flow  was  underturned  over  70  percent  of  the  span  and  the 
maximum  underturning  was  about  five  degrees. 

Since  the  metal  trailing  edge  angle  of  the  IGV  could  be  adjusted  readily  in 
Q3DFLO'8I,  attempts  were  made  to  simulate  the  measured  IGV  exit  flow  angle  in  Case 
2  by  adjusting  the  IGV  metal  angle.  This  case  is  denoted  as 
•  Case  5:  Case  2  with  Simulated  Measured  IGV  Exit  Flow  Angle 

Cases  I  and  2  are  shown  replotted  in  Figure  25  with  the  prediction  for  Case  5.  Case 
5  shows  a  fairly  constant  improvement  over  the  measured  range  of  flow  coeflicients.  It 
reduces  the  difference  between  prediction  and  measurement  of  the  pressure  rise 
coeflicient  by  about  40  percent  at  the  on-design  flow  coeflicient  but  has  progressively 
less  effect  for  higher  flow  coeflicients. 

J.  OBSERVATIONS 

From  the  above  studies,  it  is  clear  that  it  is  desirable  to  have  provision  for 
introducing  tangential  aerodynamic  blockage  as  input  data.  The  presence  of  the  struts 
ahead  of  the  IGV'  is  typical  for  an  axial  flow  compressor.  Its  effect  on  the  flow  velocity 
at  the  inlet  of  the  IGV  must  be  considered,  and  this  can  be  accounted  for  by  using  the 
tangential  aerodynamic  blockage  factor. 

The  prediction  of  the  performance  near  the  design  condition  is  much  better  than  at 
ofl-design  conditions.  An  accurate  description  of  the  flow  conditions  at  the  inlet  of  the 
compressor  is  essential  if  the  code  is  to  give  good  predictions  at  following  stations. 

The  prediction  of  the  IGV'  exit  flow  angles  by  Q3D FLO'S  1  differed  from  the 
measured  values.  As  a  result,  the  predicted  flow  downstream  of  the  IGV  exit  was  not 
representative  of  the  physical  flow.  It  is  noted  that  the  code  does  not  contain  a  separate 
correlation  for  deviation  angle  from  IGV's,  but  treats  the  IGV's  as  a  stator  (compressor) 
blade  row. 
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V.  THE  INVISCID  BLADE-TO-BLADE  SOLUTION 


A.  BACKGROUND 

Q3DFL0'81  has  a  separate  in  viscid  code  to  solve  for  the  velocity  and  pressure  fields 
on  the  axisymnietric  streamsurface  between  two  adjacent  blades.  The  blade-to-blade 
code  requires  the  inlet  and  outlet  flow  conditions  to  be  specified.  The  code  can  be  run 
independently  by  entering  the  inlet  and  outlet  conditions  manually.  Alternatively, 
Q3DFLO'81  can  optionally  couple  the  throughflow  code  with  the  blade-to-blade  code. 
In  this  case,  Q3DrLO'8I  would  first  compute  the  results  for  the  throughflow.  It  then 
transfers  the  necessary  data  internally  to  the  blade-to-blade  code,  for  the  specified  blade 
row,  and  continues  with  the  blade-to-blade  computation. 

The  blade-to-blade  code  was  applied  to  Cases  2  and  5  (of  the  preceding  chapter)  for 
flow  cceflicients  of  0.61  and  0.67  and  the  inviscid  blade-to-blade  solution  for  the  second 
rotor  at  the  tip  section  was  generated. 

B.  INPUT  DATA  AND  SOLUTION 

The  blade  profile  and  the  computational  mesh  for  the  rotor  tip  calculation  are 
shown  in  Figure  26  and  Figure  27  respectively.  A  total  of  25  stations  was  used. 
Stations  6  and  20  are  the  inlet  and  outlet  stations  of  the  rotor  blade  passage.  There  are 
two  options  for  specifying  the  outlet  flow  condition.  Either  the  Kutta  condition  is 
imposed  at  the  trailing  edge  or  the  e.xit  relative  flow  angle  must  be  specified.  The 
solution  presented  herein  is  based  on  the  use  of  the  Kutta  condition. 

As  an  e.xample  of  the  blade-to-blade  computation,  the  results  are  presented  here  for 
Case  5  at  a  flow  coeflicient  of  0.61.  The  computed  iso-pressure  lines  on  the 
a.xisyrmnetric  stream  surface  at  the  tip  and  the  pressure  coefficients  on  the  suction  and 
pressure  sides  of  the  blade,  are  shown  in  Figure  28  and  Figure  29  respectively. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  finite  element  code  Q3DFL0'81  was  applied  successfully  to  two  different  axial 
compressor  geometries.  First,  application  to  a  oenchmark'  compressor,  for  which  very 
complete  flow  survey  data  were  available,  showed  the  importance  of  accounting  for 
peripherally  non-uniform  flow  (blockage)  effects.  Second,  application  to  the  NFS  axial 
research  compressor  demonstrated  clearly  the  difficulties  that  are  inherent  in  comparing 
code  predictions  with  experimental  data  when  the  experimental  information  is  not  fully 
comprehensive,  and  when  the  measured  flow  is  not  valid  as  an  inlet  boundary  condition 
for  the  code. 

In  spite  of  the  complex  nature  of  the  flow  within  the  NPS  compressor  test  rig,  when 
a  valid  inlet  boundary  condition  was  derived  from  the  measured  data,  the  throughflow 
code  of  Q3DFL0'81  predicted  the  pressure  rise  fairly  well  near  the  design  operating 
condition.  An  analysis  of  the  measurements,  in  order  to  derive  a  proper  boundarx’ 
condition  for  the  code,  was  required  in  order  to  achieve  agreement.  At  olf-design 
conditions,  the  code  was  not  as  successful.  For  flow  near  to  the  design  condition,  the 
inviscid  blade-to-blade  code  was  used  successfully  to  generate  the  iso-pressure  lines  and 
the  blade  surface  pressure  distributions  for  the  tip  section  of  the  second  rotor. 

From  the  abo\e  experience,  the  following  were  found  to  constrain  the  application 
ofQ3DFLO  81: 

•  The  code  does  not  accept  a  tangential  aerodynamic  blockage  factor  as  an  input  for 
the  throughflow  prediction.  This  did  not  allow  the  strut  wakes  at  the  IGV  inlet 
measurement  plane  to  be  introduced  to  the  throughflow  in  the  manner  that  they 
occurred  physically. 

•  The  code  allots  the  cndwall  blockage  factor  at  the  inlet  station  equally  between  the 
hub  and  the  case  walls.  This  was  not  representative  of  the  flow  in  the  compressor 
test  rig.  The  blockage  at  the  inlet  was  largely  due  to  the  boundary  layer 
displacement  thickness  at  the  case  wall.  The  hub  had  a  negligibly  thin  boundary 
layer. 

•  The  code  only  permits  the  meriodional  mesh  lines  to  be  clustered  either  towards  the 
case  wall  or  the  hub.  In  applying  Q3DFLO'81  to  the  benchmark  data  set  and  the 
compressor  test  rig.  it  was  found  that  it  would  be  desirable  if  the  meriodional  mesh 
lines  could  be  clustered  towards  both  the  case  wall  and  the  hub.  This  is  because 
abrupt  changes  in  the  flow  conditions  tend  to  occur  at  both  the  case  wall  and  the 
hub. 

•  The  code  limits  the  capacity  of  the  mesh  to  600  nodes  and  the  input  values  for  each 
variable  such  as  the  total  pressure  at  the  inlet  station  to  10  points.  A  larger 


capacity  is  required  to  model  a  multi-stage  geometry  using  grid  lines  clustered  near 
the  walls. 

•  The  use  of  two  dimensional  cascade  data  for  calculating  losses  across  the  rotor  and 
stator  blade  rows  was  shown  not  to  be  accurate  for  a  blade  row  with  corner  stall. 
It  is  therefore  advisable  to  establish  whether  stall  regions  are  present  in  any  new 
compressor  to  which  the  code  is  applied. 

•  The  code  did  not  predict  the  exit  flow  angle  of  the  inlet  guide  vane  correctly. 
Although  this  was  overcome  by  adjusting  the  geometrv’  of  the  inlet  guide  vane  in 
the  code,  the  correction  required  that  the  actual  flow  angle  be  known.  This 
information  would  not  generally  be  available  in  a  normal  application  of  the  code. 

The  use  of  the  code  in  the  tip-clearance  investigation  on  the  NPS  compressor  is  seen 
to  be  limited  by  the  above  constraints,  but  also  by  constraints  on  the  experiment  itself 
Tor  example,  the  spanwise  distribution  of  the  strut  wake  profile  at  the  inlet  measurement 
plane  needs  to  he  measured  accurately.  Although  the  assumptions  of  uniform  wake 
thickness  and  wake  diffusion  on  the  constant  radii  axisymmetric  stream  surface  worked 
fairly  well,  an  accurate  representation  of  the  inlet  velocity  near  the  case  wall  is  required 
to  obtain  an  accurate  prediction  of  the  pressure  rise,  and  a  inore  precise  inviscid  solution 
at  the  rotor  tip. 

Unfortunately,  the  compressor  is  not  in  general,  equipped  for  peripheral  flow 
surveys,  and  the  required  modification  to  the  heavy  case  wall  would  be  difficult  and 
expensive  to  make.  It  is  clear  however,  that  peripheral  survey  data  and  the  means  to 
input  such  data  into  the  code,  are  essentia). 

Finally,  it  is  noted  that  the  present  assessment  was  made  using  an  early  edition  of 
the  code.  Later  versions  of  the  code  might  be  less  constrained. 
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Figure  1.  Stream  Siirt'aces  Tliroiigli  an  Annular  Blade  Run 

Source:  Varva  M.  11.,  Aeroilierniodynatnics  and  Flow  in  Tiirboniachines,  Wiley,  New 
■^'ork,  .May  1974. 
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Figure  3.  Complex  Flow  Phenomena  in  a  Compressor  Rotor 
Source:  Naval  Postgraduate  School  Report  NPS67-84-022,  Computational  Models  for 
Turbomachinery  Flcfvs,  by  Hirsch,  C.,  Dec  1984. 


Schematic  Meriodioiial  Section  of  a  Compressor  Stage 
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Figure  5.  Geometry  of  the  Benchmark  Compressor 


Figure  6.  Computational  Mesh  fi 
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Figure  10.  Spamvise  Distribution  of  the  Total  Pressure  Coefficient  at  Stations  7  and  13 


Figure  11.  Spanwise  Distribution  of  the  Relative  Total  Pressure  Coefficient  at 
Stations  1  and  7 


Figure  12.  Spamvise  Distribution  of  the  Static  Pressure  Coefficient  at  Stations  7 
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Figure  13.  Spamvise  Distribution  of  the  Axial  Velocity  Component  at  Stations  7  and  13 


Figure  14.  Spanwise  Distribution  of  the  Absolute  Flow  Velocity  at  Stations  7  and  13 
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Figure  15.  Spanwise  Distribution  of  the  Relative  Flow  Velocity  at  Stations  7  and  13 
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Figure  16.  Spanwise  Distribution  of  the  Relative  Flow  Yaw  Angle  at  Stations  1  and  13 
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Figure  17.  Spanwise  Distribution  of  the  Absolute  Flow  Yaw  Angi 


iwise  Distribution  of  the  Flow  Properties  Measured  Ahead  of  the 
Guide  Vane 


Figure  21.  Compressor  Strut  Section 
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Spanwise  Distribution  of  the  MLxed-out  Flow  Properties  Calculated 
Ahead  of  the  Inlet  Guide  Vane 
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Figure  24.  Spaiiwise  Distribution  of  the  Flow  Angle  from  the  Inlet  Guide  Vane 
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Figure  25.  Effect  of  the  Inlet  Guide  Vane  Flow  Underturning  on  the  Pressure  Rise 
Coefficient 
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Figure  26.  Profile  of  the  Rotor  Tip  Section 
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ure  29.  Pressure  Coefrkieuts  from  the  Dlade-to- Blade  Solution 


APPENDIX  A.  DERIVATION  OF  THE  RELATIVE  THROUGHFLOW 

EQUATION  FOR  Q3DFLO'81 


The  steady,  relative  equation  of  motion  (Eq.  18),  is 

lEx(Vx  F)=  7V^S- +  ^  {18} 

The  two  force  terms  and  Ff  are  eliminated  from  Eq.  (18)  by  projecting  Eq.  (18)  in  the 
direction  of  F^  x  IV  .  Since  is  opposite  in  direction  to  IK 

(f,x  iT').T^=0  {20} 

and,  clearly 


(F,xlf).F,  =  0 


{21} 


Hence,  the  balance  ofEq.  (18)  becomes 

-  (F,  X  in  .  IK  X  (V  X  K)  =  (fj,  X  IK)  niVjiS  -  V^f) 


Next,  Eq.  (22)  is  expressed  in  cylindrical  coordinates.  This  is  done  term  by  term. 

Step  1:  Expanding  /•},  x  IK  gives 

Ffy  X  11  =  2  —  F +  igiFfy^ll  r  —  Fj^^in)  +  i^{F ^i-ll  e  —  F^^qU  r)  {23} 

Step  2;  TV f^S  —  V^/  on  the  RHS  of  Eq.  (22)  becomes 

Since  the  model  assumes  axisymmetric  flow  by  averaging  the  flow  properties  across  the 
blade  passage 


1 

r 


{25} 
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and  Eq.  (24)  is  reduced  lo 


TU  Q  r?  r  •  /  V-  cS  ul  .  .  oS  cl  . 

T^rS  -VRi  =  hit  — - +  - ;— ) 

cr  or  cz  cz 

Step  3:  Using  Eq.  (23)  and  Eq.  (26),  the  RHS  of  Eq.  (22)  becomes 

(4  X  l»0  •  (7V,S  -  V,D  =  [F^W,  -  F„Wi,){T^  -  ^ ) 

Ur  Ur 

C  Z  C 

Step  4:  Since 

_  -  -  I  ,  f',  -’i;  ey, . ,  1  ,  «'■'«) 

^xi“'.T(75 - 

using  Eq.  (25), 


{26} 


{27} 


V  X  (■  = 


- 1 
h  y 


CZ 


dz 

dr  ’ 

cV, 

cr 

dr 

cr 


^irV.) 


ill 

cB 


)  {2S} 


cr 


{29} 


therefore 


-  -  Iff)  c(rl'f,)  cl'  cK 

IV  X  (V  X  I  •)  =  /,(  -  IV,( 

-  IV,  c{rV^)  U\  clrVg)  ^ 

+  ‘ei  ~  “  - - ) 


cz 


cr 


-  cVf  c»'  IVff  c(rVf)) 


cr 


cz 


{30} 


Since  there  is  no  difference  in  magnitude  and  direction  between  V^  and  IJ\  and  I’,  and 
IV^,  Eq.  (30)  can  also  be  expressed  as 

-  -  IV„  CirV's)  cll\  civ, 

(r  X  (V  X  I  •)  =  /,(  - IV,(  )) 


cr  cz  cr 

IV,  c(rV,)  IV,  c(rV,) 


+  ‘e(~~r 


cz 


>'  dr 


■) 


(31} 


"V  hi  ^  r(  )  y  2_  ) 


CZ  or 


CZ 


57 


Step  5:  L'sing  Eq.  (23)  and  Eq.  (31),  the  LHS  of  Eq.  (22)  becomes 


-  c(rV^)  clVr  clV, 

(F,  X  in  .  JV  X  (V  X  V)  =  (F,,  W,  -  F,,lV,)i  -jf-  -  IV^(  - ^  )) 

or  c i  Of 

H',  c(rl'f)  n;  c(r\'„) 

+  - (32) 


cz 


cr 


div,  din 

+  (fi, »'» -  fw  IKX ^ )  — F  -tF-  ) 


cr 


02 


Re-arranging, 


(F,  X  fi>) .  if  X  (V  X  n  =  ( ^  - -^  )(f,,  ir,' -  E,,  IE,  If',  +  f,,  H;' -  f,, 

Or  0  ^ 

+  +  fbr^^'AK) 

I  c(rl'fl)  , 

-  F,Jl\lV,  +  F,,\vl  -  F,,ll\\Ve) 

cz 

=  (  ^  -  ^  r  +  '’^V) ) 


{33} 


4-' 


mr^elVe  +  hr^^'r)  ~  f,.W9  +  Jf'r)) 


or 


cz 


{F,,{W]FlVl)-\\\{F,,lV,  +  F,,lin) 


Substituting  IV^  =  llj  +  11}  +  into  Eq.  (33)  yields 


(f,  X  n>  ir X  (V  X  T)  =  ( 4!^ - )(F,,iy^  -  iv,(f,,u,  +  F.^ir,  +  E,,n;)) 


-f 


cr  cz 
1  clrin 


cr 


-  fK(r,rin + +  F^^ye)) 

o  ^ 

civ,  eiv,  2  -  - 

=  ^ “fT  ■  ^ ~ 

I  d{rVa)  ^  -  5 

1  G{rVa)  ,  ^  ^ 

+  ±^jf-iFi„lV^-lV,{F,.in) 


{34} 
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Since  =  0  in  Eq.  (34),  the  LHS  of  Eq.  (22)  becomes 


{Ft,xW).  W  X  (V  X  F)  = 
,  cW,  clV, 


1  ,  1  c(ry,) 

r  or  ^  ''  F,,  dz 


Step  6:  From  Eq.  (27)  and  Eq.  (35),  Eq.  (22)  becomes 

.  F,z  1  cirV,)  hr  1  5{rV,) 

cr  cz  ^  F,,  >■  dr  ^  F^  r  cz  ^ 

F  ,  F 

J  ^  bO  ^b6 

Step  7;  The  body  force  components  from  Eq.  (36)  can  be  eliminated  by  expressing  them 
in  terms  of  the  blade  surface  angles.  Figure  30  shows  that, 

-^  =  -tanj9  {37} 


Substituting  Eq.  (37)  and  Eq.  (38)  into  Eq.  (36)  and  re-arranging, 

c\V,  dn\ 

or  cz 

( U'  +  U  fl  tan  fJ)(  4^  -  7-4^  )  -f  ( IF.  tan  ?7  -  IV, ){ -  f4^  ) 
cr  cr  cz  cz 


-  —  ( tan  /)  — - 1-  tan  tj  — :: - ) 

'  or  cz 


Step  8:  Lastly,  the  streamilinction  (ip)  defined  through  the  equations 


dtp 

2n  — —  =  p2nrbW 
cr 
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and 


=  ~p2nrbW, 


{41} 


where  b  is  the  tangential  blockage  ratio,  is  introduced.  The  LHS  of  Eq.  (39)  can  be 
expressed  in  terms  of  \J/  by  differentiating  Eq.  (40)  and  Eq.  (41),  which  yields 


dr  dr  ^  prb  dr 


and 


cU\ 


cz 


d  (  I  . 

dz  ^  prb  cz  ’ 


{42} 


{43} 


Finally,  Eq.  (19),  which  is  the  radial  equilibrium  equation  for  the  rotor  flow,  is  obtained 
by  substituting  Eq.  (42)  and  Eq.  (43)  into  Eq.  (39),  yielding 


?  .  1  cyji  e  \  c>Jj 

cr  prb  cr  cz  prb  cz 


)  = 


(fr'^+  JVfl  tan/?)( 


cl 


-  T 


ds 


cl 


2  ‘  +  (11  tan  —  fF;.)( 


cz 


-  -y  ( tan  /? 


c(rl\ 


cr 


+  tan  rj 


IF' 

c(rr,) 

dz 


) 


{44} 


60 


X 


61 


Figure  30.  Relationship  between  Blade  Surface  Angles  and  Body  Forces 


APPENDIX  B.  CALCULATION  OF  THE  MIXED-OUT  CONDITION  FOR 

THE  STRUT  WAKES 

The  process  of  the  diffusion  and  mixing  of  the  strut  wake  on  an  axisymmetric  stream 
surface  at  a  constant  radial  position  can  be  modelled  as  constant  area  two-dimensional 
mixing  of  two  streams  initially  separated  by  the  wake  thickness  (d'‘)  as  shown  in 
Figure  31.  Since  the  flow  Mach  number  is  less  than  0.2,  the  velocity  and  total  pressure 
after  mixing  can  be  determined  by  assuming  the  flow  to  be  incompressible.  By 
conservation  of  mass, 

f){s-  d'')V\^  psV-i  {45} 

so  that 

F,  =  (1-4^)F,  (46) 

Using  the  definition  of  blockage  factor  (k) 

A=(l--^)  147} 

Eq.  (46)  becomes 

=  {48} 

By  conservation  of  momentum, 

(P,  -  P2)s  +  /H.^  -  -  psV'l  =  0  {49} 

Using  Eq.  (46),  Eq.  (49)  becomes 

(F,  -  P2)s  +  p(s  -  d^)V]  -  ps(l  -  )'k?  =  0  {50} 

Re-arranging, 

(F2-F,)  =  P(l-4-)^'t(l-(l-^))  (51} 
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Using  Eq.  (h1),  Eq.  (51)  becomes 


{52} 


From  the  definition  of  total  pressure 


p„-p,+j-pvi 


{53} 


and 


Pa  =  P2+-TP^'i 


{54} 


Using  Eq.  (48).  Eq.  (54)  becomes 


{55} 


Subtracting  Eq.  (5?)  from  Eq.  (55) 

(^^2-  =  +  1) 


Using  Eq.  (52),  Eq.  (56)  becomes 

(Pa  -  Pn)  =  -  (1  -  k))  +  \pV\[k-  ~  1) 

Re-arranging, 

ra  =  P-^-\pV]{\-2k  +  k^) 
or 

P,2-Pn-\p'']i\-kf 

Mm 


(57} 


{58} 


{59} 


Mence.  the  velocity  and  total  pressure  after  diffusion  and  mixing  of  the  strut  wakes  can 
be  calculated  by  using  Eq.  (46)  and  Eq.  (59)  respectively. 
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Schematic  for  Constant  Area  Mixing  Calculation  for  the  Strut  Wakes 
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